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Abstract. 

We present polarimetric observations of the afterglow of gamma-ray burst (GRB) 021004, obtained with the Nordic Optical 
Telescope (NOT) and the Very Large Telescope (VLT) between 8 and 17 hours after the burst. Comparison among the obser- 
vations shows a 45 degree change in the position angle from 9 hours after the burst to 16 hours after the burst, and comparison 
with published data from later epochs even shows a 90 degree change between 9 and 89 hours after the burst. The degree 
of linear polarization shows a marginal change, but is also consistent with being constant in time. In the context of currently 
available models for changes in the polarization of GRBs, a homogeneous jet with an early break time of % ~ 1 day provides a 
good explanation of our data. The break time is a factor 2 to 6 earlier than has been found from the analysis of the optical light 
curve. The change in the position angle of the polarization rules out a structured jet model for the GRB. 

Key words, gamma rays: bursts - polarization - radiation mechanisms: non-thermal 

1. Introduction measurement s and o btained upper limits llHiorthet alj Tl999: 

Co vino et al. 2002d) show that the level of polarization 
The generally accepted source of gamma-ray burst (GRB) „ „ , . , „ 

° J * ° . is generally small, presumably because the intrinsically 

afterglow emission is synchrotron radiation, produced , . , , . .. , . . t 

b J r high polarizatio n is averaged out to the few percent 

when the initial relativistic blast wa ve hits the cir- , , Tp, : T~TT T — [TTwl : ] linnnl 

. j— — ; — I observed ( Gruzinov & Waxman 1999; Gruzinov 1999; 

cumburst mailer and starts radialm.fi (Recs & Meszaros |jy[ e j ve j ev & LoebTl999T) 

If the outflow of the blast wave is collimated into a jet 
radiation is highly polarized, up to 75% (|Rybicki & Lightmatj jRhoadsll 19971 1 19991 ISari et aljfl999l) . several models predict 
ll979j), and polarization has indeed been measured for 6 changes in the degree of linear polarization from a few up to 



19921 IPaczvnski & Rhoadl 


1993; Meszaros & Reed 


1997; 


WiiersetalJ Il997l IWiiers 


& Galamal 1 1999b . Svnchrotron 



GRB afterglows llWiiers et all 1 19991 < ICovino et all Il999| 3Q% Is a rli9 99: Ghisellini & Lazzati 1999; Rossi et al. 2002). 
Rol et all 12000: Biornsson et alJl2002l ICovino et alJl2002alel Sofar, only hints for these variations have been seen llRol et all 



a^ M .j^^pay aiMM ^^^j iaa j^M aK-M M^. i i ■ w JUiOl, \JLliy lllllLS 1U1 LllCdC V til 1 ill IK H 1 > liaVC UCC11 \ 1\U1 CL dlj 

l2003al jMasetti et all 120031). See also the re- |20QQh . mainly because of the low polarization values that are 

views by i|Bjornsson||2003]) and (|Covino et alj|2003a|). These mea sured and the difficulties involved in obtaining a time se- 

Send offprint requests to: E. Rol (eveft@Ktoice.Hva.iiO ries of accurate P olarization measurements of GRB afterglows 

* Based on observations made with the Nordic Optical Telescope; One exception is possibly GRB 020405, for which|Bersier et alJ 

based on observations collected at the European Southern J2003<J) find a large variation in the deg ree of linear polariza - 

Observatory, Chile, by GRACE (Gamma-Ray Burst Afterglow tion within a short time interval (see also lCovino et all2 003b). 

Collaboration at ESO), under programme 70.D-0523(A) which cannot be reconciled with any current model. 
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GRB 021004 was localized with the wide-field X- 
ray Monitor (WXM) on board the High-Energy Transient 
Explorer-II (HETE-II) with an initial positional error of 10'. 
The position was immediately issued to the community, 
which allowed the rapid discov ery of its afterglow with the 
Oschin/NEAT robotic telescope jFoxl2002l) . 

The afterglow light curve is well covered and shows some 
deviations from a standard power-law decay, for which various 
explanations have been offered, such as variations in the burst 
energy or variations in the d ensity of th e surrounding medium 
(see for example [Lazzat i et alJ l2002l |H evl & Perna 2003; 
iNakar et alJ l2003t iDado et alJ 120031) . iHolland et all ll2003h 
measure a break in the light curve between 3.5 and 7 days after 
the burst. The redshift of the afterglow pl us host galaxy wa s 
determined to be z = 2.33 (see for example Mai ler et all20 02). 
while the spectrum consists of a complex of absorption systems 
JSalamanca et al.l2002UMirabal et alJ2OO2tlM0ller et all2002l) 
at the redshift of the host. 

Polarizati on measurement s were obta i ned by various 
group s dCovino et al.1 l2002bt IWang et all 120031 iRol et alJ 
l2002h . Here, we report on early polarimetric observations ob- 
tained by our group with the Nordic Optical Telescope (NOT) 
at the Observatorio del Roque de los Muchachos on the Canary 
Islands, and later with the Very Large Telescope (VLT) at the 
European Southern Observatory (ESO) in Chile. 



2. Data reduction and analysis 

Polarimetric observations of the afterglow of GRB 021004 at 
the NOT were performed from October 4.859 UT until 4.908 
UT (~8 to 10 hours after the burst) with the Andalucia Faint 
Object Spectrograph and Camera (ALFOSC), using two cal- 
cite plates at different orientations and a Bessel R filter. The 
calcite plate yields two overlapping images of the field-of-view 
(FOV) of the telescope, separated by about 15 arcsec. One im- 
age allows for the measurement of the ordinary ray, the other 
image gives the extra-ordinary ray, which together allow for the 
measurement of one of the linear Stokes parameters. For each 
observation, the orientation of the calcite plate was either 45 
degrees or 90 degrees. We obtained 3 pairs of exposures, each 
exposure with an integration time of 600 seconds, which al- 
lows the determination of Stokes Q (0/90 degrees polarization 
orientation) and U (45/135 degrees) for each pair. 

The polarimetry observations with the FOcal Reducer/low 
dispersion Spectrograph (FORS 1) at the VLT Antu were per- 
formed a few hours later, from October 5.151 UT until 5.196 
UT. We used a Wollaston prism with a rotatable half-wave plate 
at four different angles. The images do not overlap but are sep- 
arated by a mask covering half the FOV. Each angle allows 
a measurement of both the ordinary and extra-ordinary ray, 
which in turn allows the determination of the Stokes Q and 
U parameters. The broad-band filter applied here was Bessel 
V. The observations consisted of three sets of four exposures, 
with an integration time of 120 seconds for each of the first four 
exposures, and an integration time of 300 seconds for each ex- 
posure in the last two sets. 



All data were reduced using the IRAF 1 software suite. 
The images were first bias subtracted. The flatfielding of the 
ALFOSC images was performed with the orientation of the cal- 
cite plate for the flatfield identical to that of the target image. 
The FORS 1 images were flatfielded using flatfields without the 
Wollaston prism and half-wave plate. Artefacts introduced by 
the prism or half-wave plate can be corrected for using two ob- 
servations, one with the half-wave plate oriented at either or 
22.5 degrees, and one with the half-wave plate oriented at 45 
degrees difference. 

The ALFOSC images show filter scratches which could not 
be completely corrected for by flatfielding. To allow for the 
detection of errors introduced by such artefacts, the position of 
the afterglow on the CCD was slightly offset from the centre on 
the third set of observations, while it was centered for the first 
two sets. Observations of standard stars verified that a source 
positioned in the centre gives the correct polarization. 

To measure the flux, we used aperture photometry on both 
the ALFOSC and FORS 1 images, where the apertures sizes 
were adapted to the measured seeing. Due to the small FOV of 
ALFOSC in polarimetric mode, only a few stars could be mea- 
sured, and it was not possible to derive an accurate point-spread 
function (PSF) for the two images and perform PSF photome- 
try. 

For the determination of the Stokes parameters, only the 
relative flux of the ordinary and extra-ordinary ray is required. 
In the ALFOSC case, however, the calcite plate projects im- 
ages for the ordinary and extra-ordinary ray differently, which 
results in the PSFs being different for the two images, and, for 
small apertures, the incorrect ratio of light to be measured. To 
obtain the most accurate as well as precise result, we chose 
the smallest aperture size for which the resultant Stokes pa- 
rameters are still consistent with those measured using larger 
apertures. An aperture with a radius equal to ~ 10 pixels, or 
1 .5 times the seeing, provided the best results. This aperture i s 
somewhat larger than the one used bv lBiornsson et all (|2002), 
who performed polarimetric observations of GRB 010222 with 
ALFOSC, but their observations were performed with better 
seeing (~ 079 compared to ~ l'/2 for our observations). 

Calibration and verification of the procedures was done 
with both zero polarization standard stars, to check for any 
instrumental polarization, and high polarization stars. For 
ALFOSC we used BD +28°421 1 and BD +32°3739 as zero po- 
larization stars, and HD 204827 as high polarization standard 
star. For FORS 1, BD -12°5133 was used as high polarization 
standard star. 

We performed aperture photometry on the FORS 1 data and 
verified the results by PSF photometry. We used an aperture of 
1 .5 times the seeing. As there is no difference for FORS 1 in the 
PSF between the images of the ordinary and the extra-ordinary 
ray, this resulted in modest errors. We used the other stars in 
the FORS 1 FOV to calculate the polarization introduced by in- 
terstellar matter, assuming that the net polarization of the field 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Polarimetric results 
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Fig. 1. A plot of the Stokes vectors (Q, U) and their error bars 
(1-cr of the afterglow at different epochs, and of the field stars 
as determined by the FORS 1 observation. The afterglow is an- 
notated with the time after the burst in days. The ALFOSC 
measurements are denoted with an asterisk and the FORS 1 
meas urement with a square. We have also included data from 
the Covino et al. (2002b c) (denoted with a triangle) and from 
IWang et al.l J2003) (diamond). The field stars are plotted with- 
out symbols. The dashed ellipse encloses the surface which 
contains 68% of the field stars, according to their measured 
spread. The weighted mean (Q, U) value of the field stars is 
represented by the thick cross, where the size of the cross indi- 
cates the error in the mean. We have used the spread in the field 
stars as the error in the ISM induced polarization, rather than 
the error in the weighted mean. 

stars is zero and that the largest fraction of interstellar matter is 
between these field stars and the observer. Any resultant polar- 
ization is then caused by the interstellar matter, and was found 
to be P = (0.58 + 0.33)%, 6 = (105 + 16)° Since the spread 
in the (Q, U) values of the field stars is rather large, as can be 
seen in Figure [J we used this spread as the error in the above 
values, instead of the error in the weighted mean of the field 
stars, which is much smaller. 

After calculation of the Stokes parameters, we corrected 
for polarization induced by the ISM. The degree of lin- 
ear polarization and its position angle were calculated from 
both ISM-corrected and uncorrected Stokes parameters, af- 
ter which the polarizatio n de gree was corrected following 
IWardle & Kronbergllll974h and lsimmons & Stewart! lll985h for 
bias resulting from the fact that P is a definite positive quantity. 

3. Results 

For both ALFOSC and FORS 1, we obtained 3 sets of polar- 
ization measurements. The results from FORS 1 are constant 
from one set to another, showing no variations on this short 
time scale (15 to 16 hours after the burst). We therefore used 
the summed image to obtain a higher signal to noise ratio and 
one final result. For ALFOSC, the three sets do not show ev- 
idence for such constancy, and we calculated the results sepa- 
rately (see TableHJ. Since the polarization during the ALFOSC 





ALFOSC 


ALFOSC 


ALFOSC 


FORS 1 




set 1 


set 2 


set 3 


sum 




At (days) 0.37 


0.38 


0.40 


0.67 




no ISM correction 










P(%) 1.17 ±0.46 


1.73 ±0.51 


0.15 ±0.49 


1.29 ±0.13 




<9(°) 184.2 ±11.4' 


166.4 ±8.1 


129.8 ±41.0" 


121.8 ±2.8 




with ISM correction 










P(%) 1.72 + 0.56 


2.09 ± 0.60 


< 0.59 


0.75 ± 0.42 




e(°) 187.7 ±8.3' 


173.0 ±7.9 




132.5 ± 13.9 



I We added 180 to the value of the angle for clarity. 

II The very low value for the degree of linear polarization makes the 
value for position angle very insecure, which is reflected in the 
large error. See also the text for comments on this data point. 



observation could have been variable, the assumption of con- 
stant polarization, needed to calculate P and 6 from the separate 
measurements with the calcite plate at 45 and 90 degrees, is not 
valid. To see whether the results in Table [2 are still represen- 
tative of the polarization at the times of observation, we paired 
the observations differently, obtaining two more sets. The re- 
sulting values for P and 6 are consistent with the first two sets 
of the original values, that is, they are intermediate values. The 
results for the third set are not consistent with the previous two 
results and either show a very rapid change in the polarization, 
or are due to an artefact in the data. The latter could result from 
to the aforementioned defects in the filter, caused by the after- 
glow positioned at a faulty position on the CCD (the first two 
measurements do not suffer from this, as outlined above). 

When comparing the ALFOSC data with the FORS 1 data, 
we assume there is no significant difference in the polarization 
due to the different filters (V and R) we have used. lWang et alJ 
( 2003) mention a wavelength dependent change in the polariza- 
tion for their spectropolarimetric measurements, but this only 
occurs below « 405 nm and would not affect our comparison. 

We have plotted the resultant degree of linear polarization 
and position angle as a function of t ime in Figure El where 
we have also included the results by ICovino et al.l (2002b c) 
and lWang et alJ J2003I) . From the data uncorrected for ISM po- 
larization, we see rapid changes in the degree of polarization 
between 8 and 10 hours after the burst. The polarization mea- 
sured by FORS 1 could be entirely due to the ISM polariza- 
tion, since the spre ad in the ISM polariz ation is rather large, a s 
also remarked bv ICovino et alJ J2002cl) and IRoI et alJ (|2P02). 
However, the change in the p olarization from our F ORS 1 data 
point to the one measured bv ICovino et alJ (2002c) shows that 
the polarization is at least partially intrinsic to the afterglow. 

A change in the position angle and degree of linear polar- 
ization is entirely due to the afterglow, assuming that the po- 
larization of the field stars and ISM is constant in time. The 
first two ALFOSC measurements are consistent with having a 
constant position angle; the later FORS 1 point shows a change 
by about 45° in the position angl e at a 5 sigm a level. Inclusion 
of the measurement bv lCovino et alJ J2002ch even indicates a 
change of about 90° from 9 to 89 hours after the burst. For clar- 
ity, we have plotted the Q and U Stokes parameters for all the 
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Fig. 2. A plot of the degree of linear polarization P and the 
position angle 8 as function of time. ALFOSC measurements 
are indicated with an asterisk, the FOR S 1 measurement with a 
square. We have incl uded the da t a fromlCovino et al.l (EoQ2b c, 
triangles) and from lWang et alJ (|2003l diamond). The dashed 
lines give the + l-<x ranges for the ISM polarization. 



data together with the Stokes parameters for several field stars 
in Figure [0 

To estimate the significance of this change, we calculated 
the probability that the measured values originated from one 
constant value, by calculating the y 1 values for Q and U, 

x\ = tJ j , where x — Q or U, and x is the weighted mean 

of the 7 available measurements for the corresponding Stokes 
parameter. Applying an F-test to the resultant x 1 gives the re- 
quested probability. We also applied this procedure with x be- 
ing the average ISM polarization value, and cr x the correspond- 
ing spread therein. All probabilities are small. There is a 1.2% 
chance that U can entirely be attributed to ISM polarization, 
but this probability for Q is almost zero, as is then the proba- 
bility for P. The probability that either U or Q belongs to one 
average value is less than 10~ 6 . 

4. Discussion 

Several models explain changes in the polarization angle in the 
context of jetted outflow of the gamma-ray bu rst ejecta. The 
break seen in the light cu rve of GRB 021004 jHolland et alJ 
2003; Bersier et al. 2003b) is indicative of such a jetted outflow, 
though the many bumps in the l i ght cu rve hamper the detection 
of such a break. iHolland et alJ J2003I) also exclude a spectral 
change due to the passing of the cooling frequency v c through 
the optical as a cause of this break. 

We now proceed under the assumption that the change in 
the polarization is caused by a j etted outflow . 

The model proposed by iRossi et all d2002t) explains 
changes in the polarization using a structured jet: a jet with 
brighter (and possibly faster) core surrounded by dimmer (and 
slower) wings, with a standard energy reservoir. However, this 
model does not predict a change in the position angle and is 
thus ruled out by our measurements for this burst. The model 



suggested bv lSaril i 19991) and lGhisellini & Lazzatil i 19991) pre- 
dicts a 90 degrees change in the position angle, roughly around 
the time of the break. The precise moment of this change, as 
well as the maximum observable polarization, depends on the 
ratio of 6>/6>o, where 6 is the angle between the jet axis and 
the line of sight, and 9o is the initial opening angle of the jet. 
From our measurements, we deduce that the break time, where 
9o ~ 1 /T (F being the bulk Lorentz factor), is somewhere be- 
tween 10 hours and 1 day after the burst. The dependence of the 
change i n polarizatio n o n t is via T. If T(f) is differen t than as- 
sumed in lSarilJl999h and lGhisellini & Lazzatil (fl 999). which is 
not unlikely in view of the complex behaviour of the light curve 
and its explanations, the time of the break will be different than 
estimated above. Color changes have also been seen in th e op- 
tical part of the energy distr ibution of the burst JBersier et alJ 
2003n: lMafheson et all2003l) . Such behaviour requires detailed 
models for the polarization, whereas we have here used the 
general model for a smooth afterglow behaviour. It is likely, 
however, that these color changes are of little influence to the 
models described above, since they were seen past one day, 
while the largest c hange in our data is before one day. 

lHollandetal] (2003) find a jet break time of % — 6 days af- 
ter the burst, from fitting a broken power law to the data . These 
esti mates are in st ark contrast with our findings. However, 
IHolland etai1 J2003) find that the break is gradual and occurred 
over a period of 3.5 to 7 days after the burst; they also note 
that their estimate of the break time might be too high, pos- 
sibly putting the break around 2 days. With this latter value, 
our data would agree more with the jet model for polarization. 
The various bumps in the light curve might further obscure the 
detection of an early time jet break. 

The smooth and gradual break in the light curve would 
mean that the ratio G/6q is high and we are viewing the je t 
close to the edge (see Isarill999t iGhisellini & Lazzatilll999l) . 
However, this should then give rise to significantly higher val- 
ues in the degree of linear polarization than measured, unless 
we observed close to the moment where the position angle 
changed by 90 degrees. That could then also explain the third 
set of ALFOSC observations, and would show that those ob- 
servations where taken very close to the jet break time. 

5. Conclusions 

Our polarimetric observations clearly show a change in the po- 
larization of GRB 021004, most distinct in the position angle. 
The latter changes by 45° bet ween 9 and 14 hours, a nd the in- 
clusion of a later data point bv lCovino et alJ (|2002c) indicates 
even a 90° change over a 3.5 day period. Within the currently 
proposed GRB jet models, this would mean that we are looking 
at a uniform jet, with a break time of % ~ 1 days after the burst. 
This i s in contrast with the result obtained by IHolland et alJ 
ll2003l) . who obtained th ~ 6 days, but with a large spread in this 
value (« 3.5 to 7 days), which could still be reconciled wit h our 
finding s. The structured jet model as proposed bv lRossi et alJ 
(2002) is ruled out by the fact that this model does not predict 
a change in the polarization angle. 
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